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CHAPTER  1 


BASIC  PRINCIPLES  OP  DIGITAL  BBAMPORMING 


1.0  U 


CM  A' 
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This  pipir  will  concentrate  oa  digital  beaafoneiag  la  the  reoeive 
node  for  uniformly  spaced,  liaaar  array  aataaaaa.  While  a  baaapattara 
for  a  traaaaittad  aigaal  la  daaoribad  by  ita  powar  diatribatioa  la 
tiptoe,  tha  raoaivad  baaapattara  la  aaaa  aa  tba  antenna  raapoaaa.  Tba 
two-fold  ala  of  tba  raoa iver/beaaf oraar  la  to  optiaisa  tha  daairad 
aigaal  and  ainlaixa  uadeairad  iatarfaranoa  (i.a.,  noiaa,  clutter, 
jaaaing)  [1].  A  basic  baaaforaer  and  ooatrol  aystaa  [2]  (Figura  1-1) 
coasiats  of  tha  antenna  eleaents.  tha  baaapattara  controller,  and  tha 
baaaforaar.  An  inoideat  aigaal  wavafront  iadncea  linear  phase  errors 
(delays)  aoross  tba  array  of  aatanaa  elements  (i.a.,  the  aperture). 
Tha  aystaa  removes  these  phase  shifts  froa  the  desired  direction  and 
accumulates  the  desired  signals.  The  baaapattern  controller  determines 
weighting  coefficients  to  be  applied  to  the  input  signals  of  each 
eleaent  to  aiaiaixe  phase  and  aaplitude  errors  between  eleaents,  to 
steer  the  beaa  in  the  direction  of  a  desired  signal,  to  control 
aidelobe  levels  and  null  placements,  and/or  to  produce  aultiple  beams 
for  siaultaaeous  tracking  of  several  targets.  The  beaaforaer  combines 
the  input  signals  and  the  weighting  coefficients  to  produce  the  input 
beaa  for  radar  processing. 

The  aaia  advantages  of  digital  beaaforai^g  are 


ANIFNNA 


Figure  1-1.  Basic  beamformer  and  control  system. 
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.  Aatm«  self-calibration  for  ultra-low  side  lobes, 

.  Improved  adaptive  be aaf orating, 

.  More  closely  spaced,  multiple  beams, 

.  A rray  element  pattern  correction,  and 
.  Superresolution. 

The  major  disadvantages  are  circuit  complexity,  oost  and  processing 
time . 


1.1  ANTENNA  SELF-CALIBRATION  FOR  ULTRA-LOW  SIPELQBES 

To  apply  calibration  techniques,  a  test  signal  is  applied  to  the 
input  receivers,  on  an  occasional  basis,  and  then  measured  internally 
to  offset  channel  matching  errors.  This  insures  that  the  required  beam 
shape  and  pointing  angle  (as  measured  from  the  perpendicular  axis  of 
the  trray  face)  is  maintained.  If  the  test  signal  is  applied  to  the 
antenna  elements  (by  either  a  near  or  far-field  auxiliary  antenna  or 
precise  coupling  lines  across  the  antenna  face),  then  additionally,  the 
antenna  and  feed  path  errors  are  also  offset.  For  a  test  signal 

9  9  9  9 

applied  to  each  element,  lx.,  x.,  x  ]  =  X  ,  the  response  is 

*  »  n 

9  9  9  9  9  9  9  9 

[x  , ,  x  x  ]  =  X  .  A  diagonal  matrix  operator  C  is  formed 

such  that 


i 


o 

IT 


r  =  1 , 2 ......  n 


1-1 


x 
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This  aitrix  is  applied  to  the  normal  weight  vector  W  whioh  modifies  the 
input  signals,  so  that  the  inpat  besa  shape  and  pointing  angle  are  now 
aodified  by  the  weight  vector  CW.  Problems  can  arise  from  inproper 
placenent  of  the  auxiliary  antenna.  Correct  directivity  most  be 
insured  so  that  the  array  it  illuminated  without  illuminating  nearby 
structures,  which  can  return  undesired  echoes.  Also,  the  auxiliary 
antenna  must  be  physically  placed  to  prevent  it  from  being  a  source  of 
unwanted  echo,  thus  disturbing  the  field  distribution.  With 
calibration,  a  clearer  beam  is  formed  with  sharper  nulls  and  there  is 
better  sidelobe  control.  The  system  is  not  misled  by  distortion,  and 
therefore,  responds  more  rapidly  and  smoothly  as  it  steers  the  beam. 
Also,  if  individual  elements  fail,  the  self-calibration  will  minimize 
the  effect  on  the  whole  system.  Finally,  self-calibration  can  provide 
an  on-line  maintenance  check  of  individual  elements.  If  an  element's 
performance  drops  below  a  predetermined  level,  the  operator  can  be 
notified. 


1.2  APAfTIYE.  BEAMFQ&MINfi 

The  second  advantage,  adaptive  beamforming,  may  well  be  the  best 
exploitation  of  digital  implementation.  This  technique  of  effective 
sidelobe  control  and  null  steering  utilizes  adaptive  weight  formation. 
Properly  determined  weights  applied  to  the  input  element  signals  form 
null  patterns  with  very  low  sidelobes  in  the  direction  of  unwanted 
noise  sources.  The  null  depths  are  not  affected  by  the  amplitude  and 
phase  errors  in  the  element  and  channel  paths.  Barton  [2]  points  out 


s 


that  because  of  this,  adaptive  weight  control  cottnetes  with  calibration 
techniques  in  the  job  of  nnll  steering.  This  is  apparent  since 
redneing  amplitude  and  phase  errors  by  calibration  will  not  appreciably 
enhance  adaptive  weight  control.  Therefore,  it  is  not  cost  effective 
(at  least  at  present)  to  eaploy  both  techniques  in  one  system. 

The  general  idea  of  adaptive  weight  control  is  to  determine  the 
weights  necessary  for  each  element  such  that  the  beampattern  formed  has 
nulls  in  the  direction  of  interference  (i.e.,  clutter,  jammers,  etc.). 
Thus,  with  updated  returns  the  beam  controller  discerns  between 
interference  and  the  desired  signal,  so  as  to  degrade  (or  pre-whiten) 
the  unwanted  signals  and  enhance  the  beam  response  to  the  desired 
signals. 

There  are  several  general  approaches  (or  methods)  to  adaptive 
beampattern  control.  One  of  the  earliest  is  the  Applebaam-Howells 
[3,4]  (Figure  1-2)  method  which  uses  a  closed-loop  control  system 
incorporating  an  externally  provided  steering  vector  through  auxiliary 
receiving  antenna  elements.  Sidelobe  cancellation  circuits  drive  the 
weighting  networks  toward  weight  values  that  minimize  interference 
levels  at  the  output.  This  iterative  method  was  initially  intended  for 
an  analog  system.  Closely  related  is  the  Widrow  LMS  adaptive  array 
shown  in  Figure  1-3  [3].  This  time-sampled  digital  method  uses  a 
steepest  descent  method  of  solving  a  least  mean  squares  (LMS)  error 
problem,  sometimes  called  a  conjugate-gradient  method.  It  minimizes 
the  effect  of  the  interference  at  the  beamformer  output  through  an 
iterative  scheme.  The  path  of  steepest  descent  is  based  on  some 
measurement  of  the  rate  of  change  of  output  interference  as  a  function 
of  the  in-phase  and  quadrature  components  of  the  weighting 


SUMMER 


Figure  1- 


SIGNAL 

REHJCA 


Figure  1- 


.  Applebaum-Howells  adaptive  array. 
(Taken  from  reference  4,  page  32-43) 


Widrow  LMS  adaptive  array. 

(Taken  from  reference  4,  page  32-44) 
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ooef  fioients .  A  variation  of  Widrow’s  LNS  array  \$  the  weight 
perturbation  method  proposed  by  Hudson  and  Cantoni,  which  also  obtains 
adapted  weights  using  the  steepest  descent  aethod.  Here  a  specially 
selected  perturbation  or  change  of  the  beaaforner’s  output  weights  is 
applied,  and  the  resulting  change  in  output  noise  is  measured.  This 
procedure  is  repeated  until  it  converges  to  the  weights  which  will 
ainiaize  the  effect  of  interference.  Unlike  in  the  Vidrow  and 
Applebaum-Howells  methods,  time  saaples  of  the  input  element  signals 
are  not  used.  The  biggest  problem  with  these  three  iterative  methods  is 
the  slow  rate  of  convergence  in  cases  of  unfavorable  distributions  in 
jammer  powers  and  directions. 

The  next  method  involves  the  direct  evaluation  of  the  optimum 
weights  from  a  measurement  of  the  covariance  matrix  of  the  interference 
signals  frou  the  individual  array  elements.  Compton  refers  to  this  as 
a  power  inversion  array,  while  Brennan  and  Reed  call  it  a  sampled 
matrix  inversion  (SMI)  [6]  array.  This  method  requires  digital  control 
using  complex  matrix  arithmetic.  The  input  noise  field  is  sampled  N 
times  with  the  desired  signal  absent  such  that  x^  for  j  =  1,2,  ...,n 
denotes  the  jth  sample  of  vector  x.  The  estimated  covariance  matrix  is 
defined  as 


M 


<(X(j))*(X(j)JT> 


1-2 


whioh  contains  the  values  of  interference  correlation  between  all  the 
receive  elements.  This  covariance  matrix  is  then  inverted  and  post- 
multiplied  by  a  vector  representing  the  desired  look  direction.  The 
resulting  vector  has  the  required  weights  to  minimize  interference  or 
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jammers  and  optimize  reception  of  the  desired  signals.  Probleas  with 

the  SMI  method  involve  complexity  of  the  processor  and  the  high  data 

rate  required  to  operate  in  real-time.  The  noaber  of  floating  point 

operations  (i.e.,  computer  multiplication  and  division)  is 

3 

approxiaate ty  ( 7 / 6 ) N  ,  where  N  is  the  number  of  channels.  Therefore, 
if  the  jamming  environment  changes  rapidly  with  time,  the  number  of 
elements,  N,  will  be  restricted  by  the  available  speed  of  the 
processor.  Nonetheless,  the  SMI  method  is  orders  of  magnitude  faster 
than  the  iterative  schemes  discussed  earlier.  For  large  arrays,  sub¬ 
arraying  techniques  can  be  applied  to  reduce  the  dimensionality  of  an 
NxN  array  into  N  linear  arrays,  making  a  digital  solution  feasible. 

Another  approach  to  adaptive  beamforming  is  put  forth  by  Hung  and 
Turner  [7].  Here  it  is  assumed  that  the  array  has  a  large  number  of 
elements,  N,  compared  to  the  number  of  jammers,  L,  that  the  radar  is 
designed  to  suppress.  To  suppress  L  jammers,  M  noise  samples  denoted 
by  (U1 ,  U.,  ...,  U  }  are  required,  where  L  and  M  are  approximately 
equal.  If  M  <  L,  the  first  M  strongest  jammers  are  suppressed.  The 

jammer  power  is  usually  reduced  to  a  few  dBs  above  the  white  noise 

2 

level.  This  algorithm  requires  (4M  +  6M  +  2)N  real  adds  and 

(4M2  +  8M  +  4)N  real  multiply  operations.  The  algorithm  proceeds  as 
follows : 

Step  1  Calculate  the  power  of  the  first  noise  vector 

luj2  . 

(Steps  2  to  9  use  Gram-Schmidt  decomposition  to  construct  an 


i 


orthogonal  ba.sis  set 
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{V,  >  V„ »  ...»  VM}»  M  -  numb  or  of  noise  staples.  For  the 
jaaaer  vectors,  tho  threshold  A  it  the  noist  level  below 
whioh  the  vector  it  considered  jaaaer-absent. ) 

Step  2  Fora  the  noraaliaed  bttit  vector  tt 
D,/|U  I  .  lU  |2  >  A 
1  0  *  luj2  <  A 

Step  3  Set  a  *  2. 

Step  4  Calculate  D  as 

B 


u 

m 


o  -f 

a*l 


(V  T.u  )V 
a  m  a 


Step  5  Calculate  |U  |2 

91 

Step  6  Calculate  as 

U  / lu  I  ,  lul2  >  A 

vm  -  {  B  B 

0  ,  lul2  <  A 

Step  7  If  a  =  M,  go  to  Step  10. 

Step  8  Replaoe  a  by  a  +  1. 

Step  9  Return  to  Step  4. 

Step  10  The  quiescent  weight  vector  is  represented  by  its  two 
orthogonal  coaponents  and  such  that 

is  in  the  subspace  spanned  by  the  basis  vector  V 


10 


(where  and  1°^  it  in  the  subspaoe 

orthogonal  to  V  (where  f°^V  *  0) .  Calculate  W°Q  at 

»°o  ■  ‘  L  (V.T  Vv.  • 

Stop  11  Calonlata  tha  td jotted  weight  vtotor  at 

*.  ■*V|w°ol  • 

The  Hong  and  Torner  algorithm  it  much  faster  than  previous  methods 
discussed.  For  the  example  of  a  1000-element  array  with  the  ability  to 
suppress  10  jammers,  this  method  requires  less  than  5  x  1 0 

a 

multiplications  compared  to  the  SMI's  requirement  of  10 
multiplications.  Hung's  method  is  limited  by  how  rapidly  the  jammer 
environment  is  changing  and,  therefore,  how  often  updated  weights  are 
needed,  but  appears  to  offer  the  best  approach  for  adaptive  nulling. 

The  final  adaptive  digital  beemforming  scheme  to  be  reviewed, 
based  on  the  work  of  Ward,  Hargrave  and  McWhirter,  implements  a 
systolic  array  of  parallel  processing  nodes  to  perform  their  data 
domain  algorithm  [8,S],  This  algorithm  produces  the  desired 
beampattern  without  explicitly  forming  the  adaptive  weight  vector.  If 
the  desired  signal  is  called  vector  y  and  the  weighted  array  input  is 
called  XW,  then  the  problem  becomes  the  error  minimisation  e  ■  xw-y. 
The  veotor  w  which  gives  the  smallest  vector  e  has  the  desired  weights 
(Figure  1-4).  In  the  SMI  (sampled  matrix  inversion)  method,  this 
weight  vector  is  found  by 
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forming  the  covariance  matrix,  invar ting  it.  and  multiplying  it  by  tha 
daairad  look  dirootloa  veotor.  Thia  prooadnra  it  not  oaly 
aathaaat ioally  intensive,  bat  roqairot  vary  praoisa  arithaatic  to 
pravaat  aoaarioal  iaatability  whan  tha  iapat  aat  of  liaaar  equations 
are  poorly  ooaditioaad.  Thia  iaatability  oaa  raaalt  because  tha 
eovariaaca  aatrix  iavolvaa  squaring  data  (actaally  aultiplyiag  data 
with  ita  ooajagata).  To  avoid  thia.  tha  data  doaaia  algorithm  of 
Hargrava  aad  Ward  applies  a  aariaa  of  Giveaa  rotations  to  tima  aaaplad 
data  (producing  a  Q-R  dacompoait ion) .  Tha  Givaaa  rotation  bahavaa  aa 
follows: 


»  f  » 

X  11  X12  *•*  Xln 
0  X22  *•* 

As  illustrated  iu  Figure  1-5.  a  new  row  of  data  is  applied  to  the 
triangular  array.  The  Givens  rotation  drives  the  component  on  the  left 
of  tha  pravious  data  row  to  xero.  This  modified  veotor.  reduced  in 
length,  is  passed  to  the  next  row  of  processors,  where  tha  Givana 
rotation  is  applied  again.  This  is  repeated  until  all  but  one 
component  has  been  driven  to  xero  and  the  array  has  been  effectively 
updated.  At  eaoh  step,  new  data  is  entered  at  tha  top.  The  final 
component  at  the  array  output  is  a  scaled  version  of  the  desired 
beamfonaed  output  that  would  have  bean  obtained  by  applying  the  weight 
vector  to  the  signal  vector.  Therefore,  the  weight  veotor  is  not 
derived  explicitly. 

Each  node  of  the  systolic  array  needs  to  be  an  advanced 
prograwmable  digital  signal  processor  chip  with  direot  float-point 


X11  X12  *•*  Xln 

X21  *22  •**  Sn 
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capability.  Even  to*  the  operational  bandwidth  ia  limited  (at  tha 
present)  to  10  MHi,  at  whieh  frequency  aaoh  aoda  prooaaaor  mat  pa r fora 
ita  function  within  a  tingle  Nyquist  aaapling  period  of  abont  SO  to  ISO 
at  for  real-time  oparation.  Coat  and  phyaioal  space  raatriotiona  will 
liait  tha  antanna  array  aito.  Advanoes  in  VLSI  (vary  larga  aoala 
intagration)  will  ba  naadad  for  practical  iaplanantation  of  a  ayatolio 
array. 


1.3  CLQSBLI  mm  KULTIPL6  BEAKS 

Multiple  baaaa  allow  parallel  oparation  and  higher  data  rataa  than 
possible  with  a  single  baa*,  giving  faster  assessment  of  the  desired 
region.  A  phase  array  antenna  with  N  elements  is  theoretically  able  to 
form  N  independent  beams  s imultaneonsly  from  a  single  aperture  . 
Figure  1-6  shows  simultaneous  three-beam  formation  using  three  phase 
shifters  per  element  [10],  To  form  k  beams  in  an  n  element  array  (k<n) 
would  require  nk  phase  shifters.  This  obviously  would  be  a  costly  and 
complex  system  for  other  than  just  a  few  beams.  This  antenna  is 
referred  to  as  a  post-amplification  beamforming  array.  Another  analog 
multiple  beam  antenna  is  the  Blass  beamforming  array.  Here  delay  lines 
are  tapped  at  appropriate  points  to  form  beams  at  the  desired  angles. 
One  Blass  application  used  30  miles  of  s-band  waveguide  as  delay  lines 
to  form  333  independent  beams  at  various  elevation  angles.  The  Butler 
beamforming  array  oan  form  n  beams  from  n  elements  when  n  can  be 
expressed  as  soma  power  of  two  (n  ■  2P ) .  Figure  1-7  shows  an  8- 
element.  8-beam  Butler  matrix. 
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These  three  analog  multiple-beam  formers  can  also  be  implemented 

digitally.  Once  the  iapnt  signals  at  each  element  have  been  digitised, 

the  signal-to-noise  ratio  (SNR)  is  set,  and  weights  can  be  applied  to 

form  as  many  beams  as  desired  with  no  restriction  on  their  mutual 

characteristics.  In  contrast,  RF  beamformers  must  satisfy  a  minimum 

spaoing  criterion  due  to  the  requirement  that  all  beams  be  mutually 

orthogonal  (i.e.,  ''The  average  value  over  all  angles  of  the  product  of 

one  beam  response  with  the  conjugate  of  the  other  must  be  zero'*)  [2]. 

This  arbitrary  weighting  to  form  k  simultaneous  beams  requires  4nkf^ 

real  multiplications  per  second  when  sampling  is  performed  at  f  Hz  (n 

r 

equals  the  number  of  elements  and  4  real  multiplications  are  needed  for 

each  complex  multiplication).  As  the  number  of  beams  desired,  k, 

approaches  n,  the  number  of  complex  multiplications  per  sample 
2 

approaches  n  .  Thus  for  a  large  array  the  circuit  complexity  and 

processing  time  may  be  unmanageable.  In  these  cases,  a  Fourier 

Transformation  Process  can  be  used  to  form  an  entire  set  of  n  beams 

[11,12,13,14].  If  the  discrete  Fourier  Transform  (DFT)  is  used,  which 

2 

is  analogous  to  the  Blass  beamformer,  then  4n  real  multiplications  are 
required.  A  faster,  more  economical  algorithm  called  the  Fast  Fourier 
Transform  (FFT)  reduces  the  number  of  real,  non-trivial  multiplications 
to  [2n(log2n-3)  +  8]  (for  the  Radix-2  Cooley-Tukey  FFT).  The  FFT  is 
analogous  to  the  Butler  beamformer.  There  are  other  FFT  techniques, 
such  as  the  Vinograd  FFT,  which  offer  even  greater  efficiency  under 
certain  circumstances.  The  Winograd  FFT  can  perform  the  same  transform 
using  one-fifth  the  multiplications  of  the  Radix-2  FFT.  The  Fourier 
Transform  method  of  multiple  beamforming 


number  of  elements,  N 


Figure  1-8.  Number  of  multiplications  per  range  sample  for 
different  beamforming  configurations. 

(Taken  from  reference  2,  page  272) 


F Igur  e  1  —9 .  Coupling  eflects  in  eight-element  array.  Upper  and  lower  solid  curves  show 
measured  30  dB  Chebyshev  patterns  below  and  after  mutual  coupling  correction 
Dashed  line  shows  ideal  pattern. 

(Taken  from  reference  1,  page  112) 
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differs  ia  performance  from  the  custom  beams  formed  by  arbitrary 
weights.  Arbitrary  weights  allow  arbitrary  beam  patterns  and 
directions,  whereas  the  DFT/FFT  produce  identical  beams  spaced 
according  to  the  orthogonality  constraint.  Figure  1-8  compares  the 
number  of  real  multiplications  per  range  sample  to  the  number  of 
elements  for  a  four-beam  arbitrarily  weighted  pattern,  to  (1/8)  n  beams 
using  arbitrary  weights,  and  to  the  Radix-2  FFT  with  pre-weighting. 


1.4  ARRAY  ELEMENT  PATTERN  CORRECTION 

Digital  beamforming  also  allows  the  correction  of  undesired  coupling 
effects  of  neighboring  elements  on  an  individual  element.  The  received 
signal  at  an  element  is  comprised  of  the  direct  principal  plane  wave 
and  small  signals  due  to  scattering  of  the  principal  plane  wave  on 
neighboring  elements.  For  a  linear  array,  if  the  received  signal 
vector  is  X,  and  the  desired  signal  vector  is  X  ,  then  these  can  be 
related  by  a  multiplying  matrix  C,  such  that  X  =  CX^  ,  and 

consequently,  X^  =  C  *X.  Therefore,  the  coupling  correction  can  be 
achieved  by  multiplying  tho  input  signal  vector  by  the  inverse  coupling 
matrix  C  This  results  in  improved  pattern  control  and  lower 
sidelobe  levels.  Figure  1-9  shows  the  effects  of  the  coupling 
corrector  on  an  eight-element  array,  reducing  the  sidelobe  level  from 


* 


-20  to  -30  iB. 
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1.5  SPPBRRBSOLPTION 


Resolution  is  the  ability  of  a  radar  to  distinguish  between  two  or 
aore  targets  whose  parameters  (range,  doppler,  eto.)  are  nearly  equal 
in  value.  Conventional  resolution  is  the  beaawidth  of  the  summed  bn 5®, 
which  is  limited  by  the  array  aperture  .  Superresolution  is  based  on 
the  idea  that  multiple  spatial  samples  of  the  inooming  wavefront 
provide  additional  information,  and  that  some  assumptions  abort  the 
signal  can  be  made.  Nickel  (et.al)  indicates  that  the  lower  limit  on 
superresolution  seems  to  be  two  targets  separated  by  a  quarter- 
beamwidth  [15],  He  has  data  to  indicate  that  the  limit  is  more 
determined  by  channel  mismatch  errors  than  by  pure  signal-to-noise 
ratio,  as  reported  by  Steyskal.  Superresolution  methods  involve  highly 
complex  non-linear  signal  processing  algorithms  which  will  require  a 
parallel  beimformer,  a  systolic  array  processor  for  Q-R  decomposition, 
or  a  systolic  array  processor  for  singular  value  decomposition 
(depending  on  the  actual  algorithm  selected).  Superre solut icn 
techniques  would  not  replace  conventional  resolution  techniques,  but 
provide  an  additional  refinement  when  needed. 

It*  concluding  this  section  on  the  benefits  of  digital  beamforming, 
it  should  be  stated  that  the  major  area  of  application  at  the  present 
involves  beam  control  and  sidelobe  control.  As  alluded  to  earlier,  one 
system  would  not  incorporate  all  functions  discussed. 
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1.6  PRABACTBttTSTTCS  OF  THE  INPUT  CHANNEL 

In  an  analog  system,  the  aignala  from  each  element  pass  through  a 
phase  shifter  to  remove  the  phase  delay  accumulated  across  the  array 
faoe,  and  to  steer  the  beam  in  the  desired  look  direction.  These 
signals  are  then  summed  to  produce  one  input  beam  which  is  routed  to 
the  radar  processor. 

In  a  digital  receiver,  the  input  element  signal  is  either  down- 
converted  to  some  non-zero  intermediate  frequency  (i.f.),  sampled,  and 
digitized  by  a  single  analog-to-digital  converter  (ADC),  or  down- 
converted  to  zero  i.f.  (baseband)  in-phase  and  quadrature  channels, 
sampled  and  digitized  by  two  ADC's.  Barton  points  out  that  the  extra 
hardware  in  the  second  case  is  more  than  compensated  for  by  the 
improvement  in  the  net  sampling  rate  [2].  For  a  typical  radar  pulse 
spectrum  with  a  3  dB  spectrum  width  B,  the  sampling  rate  for  the  single 
ADC  is  S.4B,  while  for  the  pair  of  ADC's  at  baseband,  the  sampling  rate 
is  1.4  B  (see  Figure  1-10).  A  diagram  of  this  input  circuit  is  shown 
in  Figure  1-11. 

1.7  CHARACTERISTICS  OF  THE  DIGITAL  BEAU  PATTERN  CONTROLLER 

Where  phase  shifters  are  used  in  an  analog  system  to  remove  phase 
errors  and  provide  beam  steering,  the  digital  system  uses  a  waveform 
digital  matched  filter  processor  to  perform  phase  rotation  arithmetic 
at  each  element.  The  digital  representation  provides  much  greater 
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0  Q7B  1  4B  frequtfnc'y 


c 


Figure  1-10.  Sampling  a  typical  radar  pulse  spectrum.  The 
pulse  in  (a)  is  shown  downconverted  to  a  non¬ 
zero  i.f.  (b),  and  baseband  frequency  (c). 
(Taken  from  reference  2,  page  268) 
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Figure  1-11.  Basic  input  circuit  with  dual  ADCs 
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flexibility  in  signal  Manipulation  and  control  due  to  algorithms  which 
oan  ba  applied. 

Two  characteristics  that  are  iaportant  in  antenna  systems  are 
dynamic  range  and  signal  bandwidth.  The  dynamic  range  of  a  digital 
beamformer  and  controller  is  determined  by  the  number  of  bits  B  in  the 
ADC.  and  the  number  of  parallel  channels  N.  Let  A  be  the  maximum 
amplitude  of  the  ADC  of  each  channel.  Then  the  total  amplitude  is  NA, 
and  the  maximum  power  is 


P  =  -iNAi2 
max  2 


1-3 


Each  channel  has  a  minimum  power. 


P 


n.min 


1-4 


and,  therefore,  the  summed  P  .  ,  assuming  the  threshold  is  set  to  where 

min 

the  average  of  half  the  channels  are  triggered  by  thermal  noisn,  is 


p  _  N  A? _ 

min  ‘  2  *  2<B-1)2 


1-5 


The  dynamic  range,  then,  is 


Pm»i  _  22*8-1) 
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-  [6(B-1)  +  l01of1QN]  dB  .  l-« 

So  it  nay  be  seen  that  the  dynamic  range  increases  6  dB  per  bit.  N 
expresses  the  gain  dne  to  "coherent  integration"  of  the  N  elements. 
Limitations  to  dynamic  range  are:  thermal  noise,  in-phase  and 
quadrature  orthogonality  errors,  quantization  errors,  sampling  time 
jitters,  and  linearity  errors. 

Signal  bandwidth  is  constrained  by  the  ADC  sampling  rate  and  the 
processing  speed  of  the  beamformer.  More  will  be  said  when  hardware 
requirements  are  exrmined. 
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CHAPTER  2 

MATHEMATICAL  BASIS  OF  BKAMFORMING 


2.0  iNTMPPCTIQli 

This  chapter  vill  develop  the  nathematlcal  basis  of  forming  a 
single  beam,  steering  the  beam,  mahing  beam  corrections,  and  forming 
multiple  beams  [16].  Signal  representation  in  the  digital  domain  vill 
also  be  examined. 


2.1  FORMING  A  SINGLE  BEAM 

Consider  beamforming  with  the  2-element  array  shewn  in  Figure  2-1. 
A  far  field  wavefront  (i.e.,  planar)  arrives  at  each  sensor  with  a 
time  delay  at  sensor  0,  relative  to  sensor  1,  of:  ?  =  dsicO/c  *= 
dsinO/Xf  (sec),  where  c  =  speed  of  light.  The  phase  shift  of  the 
voltage  at  sensor  0,  relative  to  sensor  1,  is 

0  m  2ndllB$  (radians)  .  2-1 

Let  the  received  voltagg  at  sensor  1  be  e^wt.  Then  the  normalized 
voltage  at  the  summer  output  is 

v  -  +  ejUt"d)  =  .  2-2 


% 
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To  find  M,  this  equation  117  be  manipulated  as  follows: 


,-i»  .  ,J-«  (,-J*  ,  .-!«) 


Subtract  lag  the  first  equation  froa  tbs  second  provides 


-j  t) 

e  cancels,  leaving 


t-J2»  ~  1 

-it  -  1 


v  “  e  ’ 


ritn  (e-j0/2  „  #j0/2)  * 


jut  . 
v  =  eJ  ( 


sin0/2 


)  ,  with  phase  angle  (-0/2), 


and,  finally. 


sin0/2 


This  can  be  expanded  to  N  elements  as  follows: 


jut  J (ut-0)  .  j (ut-20)  .  .  j (ut-(N-l)0) 

O  '  Q  *  0  ’  »  «  t  +  ®  1 


,J“*U  ♦  ,-J*  *  ,-JW  *  ...  . 


i9 


ti£  tift  tiC 
-j  2,  -j  2  j  2X 

e  j  (t  J  -  a-1  ) 

rje/2(rj»/2 .  #J0/2) 


*  ^sin  J/2^  with  phase  angle  (N-l)0/2  ,  and 


|v|  , 

,V|  sin  »/2 


The  aaxiena  |v|  occurs  at  0*0.  Applying  L'Hopital's  rule, 
Ivlmax  *  N.  Therefore,  the  normalized  voltage  is 


M  -  1/N 


The  beampattern  is  defined  as  the  square  of  the  normalized  voltage 


response  b(0)  versus  6,  so  that 


b(«> 


m!  .  u»2-w/a  ,  „h.r. 

K  N  siu  («/2) 


-  W-  ,U« 
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clearly  shown  in  Table  2-1. 


0  for 
b(«) 


e  for 
b<$)  , 


«  for 
b(«),Jr 


+50.9  /N 


2 

90° 

0° 

+30° 

+25.45® 

3 

42° 

0® 

+18.15® 

+16.95® 

5 

23.5° 

0° 

+10.4° 

+10.18® 

7 

16.6° 

0° 

+  7.35® 

+  7.27® 

12 

9.6° 

O'* 

+  4.25® 

+  4.24® 

50 

2.3° 

0° 

+  1.02® 

+  1.02® 

2.2  BEAM  STEERING 

If  the  signal  wavefront  arrives  at  the  array  with  an  angle  6 
greater  than  mock  the  signal  is  lost.  Therefore,  it 

becomes  necessary  to  steer  the  beam  to  receive  maximum  signal  voltage. 
This  is  accomplished  by  adding  phase  shifts  (or  weights)  to  each 
element  so  that  all  signals  reach  the  summer  in  phase,  producing  the 
maximum  input  voltage  (see  Figure  2-3). 

The  summed  voltage  is 


^  1  o  -j(i-n+l)0  ,  a  2nd  .  A 
>  S.e  J  for  0  =  sin6  . 

r=0  A 
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Now,  ■ 
Therefore, 


.JCt-tC-D-il#).  ,a4  Ci  . 
the  a  turned  voltage  of  their  products  is 


j(o»t-(a-l)0)  j(n-l)0  ,  jUt-(n-2)0)  j(n-2)0  x  jwt  jO 

V  =  0  ©  t  ©  ©  -r,*,  *  ©  ©  t 

v  *  Ne^*1**  ,  2-12 


tad  the  beampattern  respoase  is 


««  .h£ .  4 .  x . 

NZ  NZ 
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Heace,  the  beam  is  steered  to  compensate  for  the  angle  0.  and  is  now 
’'looking'*  in  the  desired  direction. 


2.3  BEAM  CORRECTION 

Just  as  the  beam  was  steered  by  applying  phase  shifting  (or 
beamsteering)  coefficients  to  the  signals  at  each  element,  corrective 
coefficients  may  be  applied  to  improve  other  aspects  of  the  received 
signal  as  well. 

A  target  in  the  noar  field  results  in  a  return  wavefront  which  is 
not  planar,  but  curved.  The  center  of  the  curved  wavefront  arrives 
earlier  than  the  rest,  causing  distortion.  The  closer  the  target,  the 
more  pronounced  is  the  error.  A  phase  lag  correction  coefficient  based 
on  target  range  can  be  applied  to  minimize  this  error.  Let 
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0^0  -  -k[<R2(t)+*i2>1/2  -  R(t)]  - 


-  2R(t)  if  xi<<R(t)* 


2-14 


where 

k  »  2n/'?.  , 

Xj  =  id  (i  *  ith  element  from  center)  , 

R(t)  =  Range  of  Target 

This  correction  is  referred  to  as  range  focusing. 

Weighting  coefficients  are  also  used  to  control  the  sidelobes 
of  the  voltage  response  of  each  steered  beam.  These  weights  adjust  the 
gains  of  the  receiving  elements  to  produce  low  sidelobes. 

The  received  signal  altered  by  range  focusing,  beamsteering,  and 
side lobe  shading  is 


B  ■■= 


WiV 


j2R(  t) 


-j(i-n+l)0 

© 
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This  is  shown  in  Figure  2-4. 


2.4  MULTIPLE  BEAM  FORMING 


For  arbitrary  weighting,  multiple  beams  are  formed  by  applying  tue 
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Figure  2-4.  Single  frequency  beamformer  output  (steered, 
focused,  and  shaped). 

(Adapted  from  reference  16,  page  9) 


B(t ,9) 
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input  signal*  to  k  sots  of  weights,  where  k  equals  the  number  of 
beams  desired,  and  then  suaaing  the  results  in  k  summers  to  produce 
k  beams  at  the  output. 

For  orthogonal  multiple  beamforming,  the  following  development 
applies. 

If  the  phase  delay  of  the  input  signal  is  represented  relative  to 
sensor  0  rather  than  sensor  n-1,  then  the  input  signal  at  any  element 
sensor  is 


=  +  id) 


2-16 


Thus,  the  beamsteering  coefficient  is  now  in  the  form 

C4  »  e"jl®  .  2-17 

The  phase  delay  0  can  be  represented  by  a  time  delay  x  by  the  following 
transfer,  action: 


d  =  sin0  ,  where  X  *  c/f  ,  so  that 

0  =  ~ n ^  sin0  *  w(^  sin0)  .  2-18 

o  c 

J  Q 

Letting  x  *  *  sin0,  where  c  =  3  x  10°  m/s,  then 
c 

-  j  iuit  _  „ . 

Ci  ~  c  =  e  J  .  2-19 


The  beaaformer  output  from  Figure  2-4  was 
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■(*.»  y  ,mt)  ■ 

T-0  1  1 


If  the  input  signal,  sidaloba  shading,  and  range  foeusing  are  combined 
to  fora  4.,  and  the  new  beaasteering  coefficient  (2-19)  is  used,  thon 


B(t,0)  *  ^  Ai« 


-j  iut 
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This  complex  summed  signal  has  the  same  form  as  the  discrete  inverse 
Fourier  Transform 


A 


where  F(k)  is  analogous  to  B(i,0),  f(n)  is  analogous  to  A^,  and 
(2n/n)nk  is  analogous  to  iwt.  This  has  importance  because  aperture 
space  can  be  transformed  to  beam  space  using  a  Fourier  transformation. 
Therefore,  an  array  with  N  aperture  elements  can  produce  N  beams.  To 
produce  N  beams  from  N  sensor  elements,  the  following  scheme  can  be 
employed: 

The  yth  beam  is  represented  as 


B(t,9Y)  -  \  A  (t)e~j(i“TY)  -  ^ 

i=0  i=0 


A.(t)e"j(2nri/N), 


y  3  0,1,2, .. . ,N-1. 
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B(t,B  /  -  AQ  +  Ax  +  A2  +  ...  +  A, 


N-l 


8(«.«1)  -  A„  ♦  *  V-J'4""1*  *  ...  *  A^'1^’ 


P(t,«N_1)  -  A 


.w^iK-JLLs  x 

+  A.e  N  '  +  ...  +  A^e  N 


0  "1 


2-23 


Etch  beta  is  centered  on  ,  where 


«T  »  sin  lNdJ  ,  X  «  0,  1,  ....  N-l. 


The  beams  formed  farther  away  from  the  boresight  axis  tend  to  have 
broader  beamwidth.  This  effect  is  approximated  by  the  relationship: 

BW,.n  *  .88d(^r)  secO  near  broadside.  2-24 

JdB  Nd 

The  endfire  beamwidth  is  larger  than  the  broadside  beamwidth  by  a 
factor  of 

Df 

"MBtE »«!»)  .  2[.88Sa/Nd)]l/J  2-2J 

3dB< Broad side) 

Beam  broadening  away  from  the  boresight  axis  increases  side  lobes  and 
limits  the  total  angle  of  spaoe  that  can  be  covered  (i.e..  the  summed 
beamwidth)  to  approximately  +60°.  Within  this  beamwidth.  it  may  be 
desired  to  have  as  many  beams  as  feasible  since  as  the  number  of  beams 
increases,  the  individual  beamwidth  decreases.  Thus,  angular 
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resolution  decreet**,  allowing  better  eeperetion  of  targets  ia  close 
proximity. 

Aa  ecoaoaioal  epplioetioa  of  the  Disorete  Fourier  Treasfora  (DFT) 
is  the  Fast  Fourier  Treasfora  (FFT).  This  is  ea  elgorltba  which  takes 
the  set  of  N  equations  (2-23)  given  above,  and  produces  N  beans, 
requiring  only  2Nt (L0GjN)-3 ] +8  real  ault ipl icat ions  (for  the 
Cooley-Tukey  Radix-2  FFT).  A  straight  forward  DFT,  in  comparison  , 
would  require  approximately  hf4  multiplications. 

It  aust  be  stressed  that  all  beamformiug  discussed  thus  far  for 
linear  arrays  is  frequency  sensitive.  The  lowest  sidelobes  occur  when 
the  inter-element  spacing  is  set  at  X/2,  in  reference  to  the  operating 
frequency.  If  a  different  frequency  is  used,  the  spacing  is  no  longer 
X/2,  and  sidelobes  increase  in  size  and  number.  The  main  lobe 
beamwidth  is  also  changed.  This  in  turn  changes  the  angular  resolution 
(for  the  worse).  Therefore,  the  phased-array,  as  described  above,  is  a 
single  frequency  system. 

2.5  QYEBALL-SISNAL  fLW 

Having  looked  at  basic  beamforming  for  an  equally  spaced  linear 
array  for  both  single-beam  and  multiple-beast,  the  signal  flow  for 
digital  beamforming,  oomplete  with  errors,  will  now  be  examined. 

The  input  signal  from  a  single  element  is 

V  (t)  *  A  cos(a»t+t  )  +  u(t)cos  ut, 
n  n  n 


2-26 
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mhere  t)  ia  the  Gaussian  noise  with  variance  a.  Aftar  conversion  to 
baseband  and  separation  into  ia-phaae  and  qaadratare  components,  tha 
aigaal  ia 


V  (t)  “  I  (t)  +  Q  (t)  ■  A  coat  +  q.  (t)  +  A  siat  qA  (t).  2-27 
a  a  a  a  a  t  n  a  u 

a  a 


Includiag  mismatches  ia  tha  ia-phaaa  aad  qaadratara  mixers,  whioh 

produce  phaaa  error  S  aad  gala  error  Z  ,  results  ia 

a  a 


I  (t)  *  A  (cost  -  1/2  S  siat  +  1/2  E  cost  )  +  q_.  (t! 
a  an  aa  an  I 


Q  (t)  -  A  (sint  -  1/2  S  cost  -  1/2  Z  siat  )  +  q_  (t).  2-29 

n  nn  an  n  n  Q 

a 


Let  these  mismatches  be  defined  as 


T*"®"  -1/2  S  sint  +  1/2  K  cost  ,  and 
A  n  n  a  a 

a 


2-30 


:  “■*  -1/2  S  cost  -  1/2  K  sint  . 
A  an  n  n 

n 


Then 


I  (t)  *  A  cost  +  +  i\T  (t)  ,  and 

a  a  a  I  *1 

n  n 


Q  (t)  *  A  sint  +  +  H_  (t) 

a  a  a  Q  Q 

a  a 


2-33 
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At  this  point,  the  signals  are  stapled  at  a  aaxiaua  rat*  o f  1.4B,  where 
B  la  the  bandwidth  or  spaetrai  width  of  tha  radar  pulse,  and  digitised. 


Lat  tha  ooaplax  waighting  tarn  ba  dafinad  aa 


W  -  IW+JQW  -  C  cos6  +  JCsinD, 


and  tha  reoeived  signal  ba 


x  *  I  +JQ  ■  S  -  AS  +  r\  , 
n  n  n  n  n  n 


where 


Sn  "  V*P<jTn)  *  ASn  "  *1  +  *  ,nd  ’’a  "  nI  +  * 

n  n  n  n 


Than 


x  *  WX  -  (I  I  -  Q  Q  )  +  j (I  Q  +IQ) 
n  n  n  w  nw  nw  wn 


Nota  that  four  raal  multiplications  and  two  additions  ata  required  for 
ona  ooaplax  ault ipl ioat ion.  Due  to  quantization  errors  in  forming 
weights,  wa  define 


W  X  -  Of  -  AW  )(S  -  AS  +  N  ) 

n  n  n  a  q  &  a 


4 
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2-37 

whin 


-IS  -  f  AS  +  W  n  -  S  AW  4-  AW  AS  -  A!  n 
a  a  a  &  a  &  a  a  aa  a  a 

-  W  S  -  f  AS  +  I*  -  SAW, 

-aa  aa  aa  aa 


and 


<|AW  |2>  ■  (1/6)  22*1  (weighting  qaaatiaatioa  error), 
a 

<|q  |2>  «  2(a2  +  (l/12)q2]  (tigael  noise  end 

quant iiation  error), 

<|AS  |2>  *  (1/4)  A2  (S2  +  I2  )  (signal  aisaatch  error), 

n  n  n  n 


2-38 


2-39 


2-40 


In  the  above  equations,  is  the  nuaber  of  bits  ia  the  weighting  word, 
and  q  is  a  single  step  of  quantized  voltage. 

The  signals  and  weights  for  the  whole  linear  array  can  be 
expressed  as  the  vectors 


S*  -  [s1.  s2 . SJ 


W 1  -  ML.  W, . W  ] 

1  *  XX 


AS  =  [AS1,  AS2,  ....  ASnl 


AW1  -  tAWj,  AW2 . AIJ 


H  m  q2 ,  .*.,  q^] 


and 


i 
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T 

X  *  ^2#  •  •  • »  • 

In  vector  for*,  the  beamformer  output  it 
Q  =  W1TX(l+p)  *  (W+AW)T(S+AS+H)<l+p), 


2  WTS  +  fTAS  +  WTH  +  AWTS  +  pWTS,  2-41 

T 

where  p  is  word  length  truncation  error.  The  first  term,  W  S, 
represents  the  desired  beamformer  output.  The  next  four  terms 
represent  first  order  errors. 

Different  weighting  distributions,  such  as  Dolph-Chebyshev 
coefficients  or  Hamming  coefficients,  may  be  used  to  control,  in 
general,  the  sidelobe  r*rpc  i.  Adaptive  weighting,  as  well,  can  steer 
beams  and  nulls,  or  reduce  sidclobes  as  a  result  of  self  calibration, 
or  modify  the  beam  shape.  Figure  2-5  shows  the  effects  of  various 
errors  in  a  digital  beamformer  for  a  Hamming  weighting  distribution  for 
N  channels. 


Thermal  noise  is  generated  by  receiver  components.  It  is  Gaussian 

2 

in  nature  and  has  power  equal  to  its  variance  or  . 

Coding  or  quantization  error  tends  to  vary  in  a  sawtooth  manner  in 
the  absence  of  thermal  noise.  As  the  level  of  thermal  noise  increases. 


maiimom  signol  l*vet 


'igure  2 


-5.  Effects  of  various  errors  in  a  digital  beamformer 
for  hamming  weighting  against  N  channels. 

(Taken  from  reference  2,  page  272) 
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the  quant izetion  error  is  reduced.  For  o>.5q,  the  error  is  negligible 
end  can  be  combined  linearly  with  the  thernal  noise.  The  quantization 
error  range  is  assumed  to  be  between  -l/2q  and  l/2q,  so  that  the  mean 
square  error  is 


[X2/q]  dx  =  <l/12)q2  . 


2-42 


Therefore,  the  net  average  noise  power  (thermal  and  quantization)  is 


<hnl2>  =  2to2  +  <l/12)q2] 


2-39 


in  each  channel. 

I  and  Q  mixer  mismatches  were  defined  in  Equations  2-30  and  2-31, 

which  lead  to  the  combined  signal  mismatch  error  power  of 

Equation  2-40,  (1/4)A2  (S2  +  K2  ) . 

n  n  n 

There  are  phase  and  amplitude  errors  between  elements,  as  well  as 
coupling  errors  (discussed  in  Chapter  One).  These  are  not  to  be 
confused  with  the  mismatches  between  I  and  Q  channels  of  a  single 
element.  Errors  may  result  if  individual  elements  malfunction  or  quit. 
It  was  described  earlier  how  a  digital  self-calibration  can  compensate 
for  these  errors.  With  adaptive  beamforming,  the  null  depths  are  not 
affected  by  these  errors. 

If  a  radar  system  is  airborne,  and  thus  on  a  moving  platform, 
linearity  errors  can  result.  Harmonic  inter-modulation  due  to  doppler 
effect  on  clutter  creates  non-linear  distortion  at  the  harmonics  of  the 
frequency  offset  error,  which  can  be  a  serious  problem  for  analog 
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systems.  Mathews  points  ont  that  digital  beamformers  are  anoh  more 
tolerant  and  that  design  criteria  may  be  relaxed  [17]. 
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CHAPTER  3 


CURRENT  TECHNOLOGY 


Systems  currently  being  developed  nsing  digital  techniques  are 
experimental  in  nature  and  are  assessing  the  directions  in  which  to 
proceed  in  order  to  reap  the  most  benefits  at  an  ' 'affordable  *  '  price. 
A  comparison  of  the  five  systems  discussed  is  found  in  Table  3-1. 


Standard  Telecommunications  Laboratories,  Limited  (situated  in 
Harlow,  Essex,  England)  jointly  with  ITT  Gilfillan  (in  JSA)  have  an 
experimental  system  for  adaptive  digital  beamforming  [2].  Their  system 
is  designed  to  reject  jamming  on  a  vertical  pattern.  A  null  is 
directed  toward  the  jammer  while  maintaining  coverage  above  and  below 
this  angle.  A  vertical,  rece ive-only,  linear  array  of  eight  horns 
covers  a  20°  azimuthal  sector  and  operates  in  the  L-band  (1-2  GHZ). 
The  rotating  test  radar  operates  normally  until  the  designated  20° 
sector  is  reached.  At  this  point,  the  no,rmal  receive  circuits  are 
switched  to  the  eight-horn  digital  receivers.  The  eight  inputs  are 
digitized,  adaptively  weighted,  and  summed  in  2  ps.  The  system 
sampling  time  is  6  ps,  providing  a  more  than  ample  margin. 


TABLE  3- 
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3.2  ELRA 

The  ELSA  experimental  phased  array  radar  is  located  at  FFM 
(Forsohungainstitut  Fur  Funk  und  Matheuatik)  in  the  Federal  Republic  of 
Germany  [18,1].  There  are  separate  arrays  for  transmitting  and 
receiving.  The  antennas  are  S-band  circular  thinned  phased  arrays. 
There  are  300  transmit  elements  with  a  1:8.2  fill  faotor  and  768 
receive  elements  with  a  1:6.25  fill  factor.  The  system  uses  analog 
phase  shifters  at  each  receive  element,  which  output  in-phase  and 
quadrature  components.  Element  signals  are  then  summed  into  48  sub¬ 
arrays  of  16  channels  each,  and  at  this  point  digitized.  The  A.D.C. 
operates  at  a  2 MHZ  rate.  Three  fixed  beams  (one  sum  beam  and  two 
difference  beams)  are  obtained.  Adaptive  weighting  and  beamforming  is 
not  done  at  present,  but  system  expansion  to  allow  for  it  is  planned. 

The  transmitter  provides  a  2-ys  tracking  pulse  and  a  10-ps  search 
pulse.  Since  the  sampling  time  is  .5  ns ,  enough  time  is  available  to 
generate  six  independent  beams  by  a  time-multiplex  method.  This  rapid 
steering  to  six  different  directions  gives  the  appearance  of  having  six 
simultaneous  beams. 

Se lf-oal ibration  is  accomplished  by  an  auxiliary  antenna  probe  in 
the  near-field.  The  phase  and  gain  of  each  element  is  evaluated  and 
the  phase  information  used  in  the  control  of  the  phase  shifters. 
Faulty  elements  with  insufficient  gain  or  ex‘cc*si?e  phase  error  are 
identified  for  future  repair.  A  complete  check  of  all  transmit  and 
reoeive  elements  takes  approximately  16  seconds. 
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3.3  GASP 

The  GEC  array-** ignal-proce sting  (GASP)  test  bed  is  an  experimental 
linear  phased  array  receiving  system  at  the  Marooni  Research  Centre  in 
Chelmsford.  Essex,  O.K.  [19].  This  system  will  have  32  elemects  then 
completed.  As  of  1985,  25  elements  have  been  interfaced,  and  the 
results  presented  by  Old  (et.al)  are  based  on  16  elements  operating  at 
S-band  (3  GHZ).  GASP  uses  a  fully  digital  implementation.  The 
received  signal  is  down-converted  to  in-phase  and  quadrature  baseband 
signals,  digitized  in  12-bit  ADC's  at  a  3-ps  sampling  rate.  The 
digital  signals  are  now  available  to  a  control  computer,  a  sample  store 
(2048  word  memory  for  off-line  processing),  and  the  real-time 
beamformer.  The  beamformer  produces  an  output  expressed  by 


Array  Output  =  ^  D  *  W  *  S  *  C  , 
Z=0  n  n  n  n 


3-1 


where 


Dq  =  Incoming  signal  from  nth  element, 

=  The  weighting  coefficient, 

S  =  Hie  beamsteering  coefficient,  and 
n 

C  3  The  correction  coefficient. 

11 

This  flexible  system  can  aocept  fixed  weighting  coefficients  (such  as 
Hanning  and  Hamming),  or  calculate  them  based  on  use r- supp  1  ie d 
parameters  (suoh  as  Taylor  or  Dolph-Chebychev) . 
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GASP  alto  incorporates  self-calibration  using  s  fsr-field  probe  or 
oonplers  pleoed  between  the  elements  end  receivers.  Phase  end  gain 
errors  ere  oelonleted  using  312  ties  samples  inpnt  to  an  FFT.  The 
results  sre  used  to  fore  end  store  the  correction  coefficients,  C  ’s. 

A 

This  highly  flexible  system  will  enable  a  wide  range  of 
experiments  in  digital  signal  processing  and  calibration  techniques. 


3.4  CQMLS-B-glfi 

This  continental  United  States  bistatic  backscatter  over-the- 
horizon  radar  is  built  by  the  General  Electric  Company  [1].  Bistatic 
radar  has  its  transmitter  and  receiver  located  in  different 
geographical  locations.  It  is  a  low-frequency  radar  (between  6-20  MIC) 
with  a  bandwidth  of  about  100  KHZ.  82  receive  elements  are  used  to 
form  four  simultaneous  beams  covering  the  transmitted  beamwidth. 
Adaptive  nulling  and  on-line  receiver  error  correction  are  performed. 


3.3  RADC 

The  Rome  Air  Development  Center  in  Bedford,  Massachusetts,  is 
testing  a  32-element  array  at  C-band  (4-6  GHZ)  [11.  This  system  has  a 
300  KHZ  bandwidth.  Self-calibration  through  feedback  to  the  element 
front-ends  should  provide  very  low  sidelobe  patterns. 
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3.6  flEMEAL  HAIDWAER  REQUIREMENTS  FOR  FULL  DTflTTAI.  TMPlJUiRNTATIQN 

The  Mia  limitations  to  digital  beamforming  art  speed.  complexity 
and  ooat.  Evan  analog  phased  array  systems  have,  for  tha  Boat  part, 
baaa  rastriotad  to  ailitary  use  whara  mission  requirements  ware  aora 
iaportant  than  ooaplaxity  and  ooat.  Much  work  ia  being  done  (such  aa 
aonolithio  microwave  integrated  oirouita  -  MNIC's)  to  reduce  the  size 
and  cost  of  these  systems,  while  increasing  component  and  processing 
speed.  As  the  overall  cost  of  phased  array  systems  is  reduced,  digital 
beamforming  will  also  become  more  attractive. 

In  the  race  Ivina  channels  of  a  digital  beamformer,  the  phase 

shifters  are  replaced  by  digital  weighting  coefficients  whioh  combine 

with  the  element  input  signal  after  analog-to-digital  conversion.  The 

ADC  and  the  number  of  elements  set  the  limits  on  dynamic  range, 

[6(B-1)+10  log  N] ,  and  the  bandwidth,  (f  /1.4).  Presently 

•iu  sampl  ing 

available  ADC's  are: 

12-bit,  5-MHZ  ADC  from  Analog  Devices  on  5''  *  T'  board  [211, 

12-bit,  2-NHZ  ADC  from  ILC  Data  Device  on  40-PIN  DIP  [221, 

8-bit.  50- MHZ  ADC  [11. 

These  would  limit  bandwidths  to  3.5  MHZ  for  12-bit  and  35  MHZ  for  8-bit 
ADC's. 

The  beam  pattern  controller  determines  the  adaptive  and  corrective 
weights  to  apply  to  the  input  signals. t  The  processing  speed  is 
determined  by  the  algorithms  used  and  the  number  of  signal  elements  in 
the  array.  An  important  consideration  here  is  how  often  the  weights 
need  updated,  whioh  for  adaptive  nulling  depends  on  the  rate  of  change 
in  position  of  airborne  jammers.  Early  digital  beamformers  will  work 
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with  ’’acceptable'',  aot  ’ ’optimal’ ’»  update  rata*.  The  choice  of 
Hung’*  and  Turner’*  algorithm  woul<*  appear  to  give  the  fastest 
throughput.  For  nulling  10  jammers  with  a  100-element  array,  the  SMI 
method  require*  a  alllion  aultiplioations.  The  Bung  and  Turner 
algorithm  requires  less  than  50,000  multiplication*.  If  updates  were 
needed  10  tiaes  a  seoond,  then  total  Multiples  are  less  than 
500,000/seo,  which  is  easily  handled  by  a  fast  dedicated  prooessor. 
Even  2,000  elements  would  need  only  10  million  operations  per  second, 
which  is  still  realizable  on  a  central  processor  [1].  Larger  arrays 
may  also  be  handled  by  sub-arraying  techniques  to  reduce  the 
computational  load. 

The  bcsmformer  must  multiply  the  individual  signals  by  the  most 
recent  weights,  steering  coefficients,  and  corrective  coefficients,  and 
then  sum  them  into  the  final  beam  output.  This  requires  (4Nkff)  real 
multiplications,  where  N  is  the  number  of  elements,  k  is  the  number  of 
simultaneous  beams,  and  f^  is  the  sampling  frequency.  One  beam  formed 
with  64  elements  at  a  1-MHZ  sampling  rate  requires 
4  x  64  z  1  x  10^  ■  256  x  10^  operations  per  second.  Present  processors 
can  handle  this,  but  not  if  thousands  of  elements  form  the  array. 
Several  alternatives  are  possible.  The  first  is  placing  a  digital 
signal  prooessor  in  eaoh  input  ohannel  to  do  all  or  most  of  the 
multiplication.  Texas  Instruments  has  for  sale  a  16-bit  x  16-bit 
parallel  aultiplier/accumulator,  the  Tact  1010-65,  which  multiplies  and 
adds  in  65  ns  [23].  This  allows  for  a  3.8  MHZ  sampling  rate  if  four 
multiplies  are  performed  at  each  channel. 

Another  alternative  is  using  a  Fast  Fourier  Transform  (FFT) 
processor,  whioh  simultaneously  forms  a  set  of  N  beams  from  N  signal 
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input s  *  requiring  {2Nt(LOG2N)-3]+8}  multiplicat ion*  pur  sample.  As 
long  it  N  <  2000.  the  FFT  takes  fewer  operations  than  forning  four 
bunas  using  arbitrary  mights.  For  thu  ossu  of  64  uluaunts.  forming  64 
bunas  requires  392  x  10^  opurntions  pur  suoond  for  n  1-KHZ  tuapling 
rntu. 

As  suun  in  Tnblu  3-1.  limits  in  prooussing  speeds  have  kept  thu 
numbur  of  army  uluaunts  low  (below  100  nftur  sub-arraying)  and  tha 
bandwidth  small  (less  than  0.5  MHZ)  on  thu  currant  experimental 
systems. 
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CHAPTER  4 


FUTURE  TRENDS 


4.0  mSIBLE.SIS.IEKS 

The  systolic  array  discussed  in  Chapter  1  nay  soon  find 
application  where  the  bandwidth  is  snail  (as  in  sonar),  dne  to 
processing  node  speeds,  and  where  the  number  of  elements  is  limited. 
For  example,  a  10-element  array  would  neud  55  signal  processors,  while 
a  50-element  array  would  need  1,465  processors.  As  progress  continues 
in  making  smaller,  faster,  and  less  expensive  signal  processors.- 
systolic  arrays  may  find  their  way  into  a  number  of  least  mean  square 
error  applications.  Some  experimental  work  is  being  done  for  adaptive 
beamforming  in  the  U.E.  [24]  and  by  RCA  in  the  U.S.  [25]. 


4.1  INIESSAXEP  email  JECBHfiLfflg 

The  push  for  single  chip  circuits  that  are  smaller,  faster,  and  of 
lower  cost  continues  in  both  silicon  and  Gallinm  Arsenide  (GaAs)  [26]. 
GaAs  logic  is  now  available  at  speeds  up  to  3  GBZ,  and  this  is 
expected  to  increese  rapidly.  HP's  C.  House  predicts  that  todays  high- 
price  6-bit  1-GEZ  ADC  will  be  a  low-cost  10-GHZ  ADC  by  1994.  NEC  will 
soon  be  offering  a  digital  signal  processor  (DSP)  chip  with  32-bit 
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floating  point  opt rations  mad  a  150  aa  oycla  tiaa.  Evan  f attar  eycla 
tiaas  and  inaxpanaiva  parallal  prooaaaiag  will  b*  naadad. 

A  aarriaga  of  aonolithic  analog  olroaita  with  high  apaad  digital 
oiroaita  aay  ba  abla  to  draw  froa  tha  baat  of  both  doaaiaa. 
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CHAPTER  5 

CONCLUSION 

Tremendous  steps  have  been  taken  in  making  digital  beamforming  a 
realizable  method  of  achi  ving  a  highly  flexible  and  responsive  system 
in  the  presence  of  rapidly  changing  noise  and  jamming  conditions. 
Although  it  may  be  the  year  2000  before  full  digital  implementation  is 
feasible  and  cost  effective,  sub-arraying  techniques  will  allow  some  of 
the  benefits  to  be  employed  before  then.  Self-calibration  may  be  useu 
sooner  to  achieve  lower  sidelobes  sinoe  it  is  not  as  restricted  by 
processing  times.  For  now,  digital  beamforming  will  be  restricted  to 
systems  of  few  elements  (probably  less  than  100)  and  small  bandwidth 
(of  about  100  KHZ). 

Three  areas  of  research  and  advancement  will  hasten  the  advent  of 
widespread  digital  beamforming:  work  on  the  speed,  size,  and  cost  of 
ICs,  the  development  of  creative  architectures  to  maximize  data 
throughput,  and  the  development  of  refined  and  new  algorithms  to  reduce 


processing  times. 
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